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Abstract 

We study the Heisenberg model in an external magnetic field on curved sur- 
faces with rotational symmetry. The Euler-Lagrange static equations, derived 
from the Hamiltonian lead to the inhomogeneous double sine-Gordon equation 
(DSG). However, if the magnetic field is coupled with the metric elements of the 
surface, and consequently, its curvature, the homogeneous DSG appears and a 
2-7r-soliton is obtained as a solution for this model. In order to obey the self-dual 
equations, surface deformations are predicted at the sector where the spins point 
in the opposite direction to the magnetic field. The model was used to partic- 
ularize the characteristic lenght of the 2-7r-soliton for three specific rotationally 
symmetric surfaces: the cylinder, the catenoid and the hyperboloid. Fractional 
27r-solitons must appear on finite surfaces, as the sphere, torus and barrels, for 
example. 
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1 Introduction and Motivation 

In the last decades, the study of relations between geometry and the physical 
properties of condensed matter systems (CMS) has attracted large attention. On the 
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one hand, the growing capacity for fabricating and manipulating nanoscale devices with 
different geometries, such as quasi two-dimensional exotic shapes like the Mobius stripe 
[I], torus [2] and asymmetric nanorings [3J, brings the possibility to develop and test 
theoretical models in several branches of CMS, e.g., nanomagnetism, nematic liquid 
crystals, graphene devices and topological insulators. On the other hand, theoretical 
models predict the appearing of collective modes, which are strongly influenced by 
the curvature of the substrate. For instance, particle-like excitations, like vortices and 
solitons, appear in several contexts, as superconductors, ferromagnetic nanoparticles, 
nematic liquid crystals, Bose-Einstein condensates [HE] and colisionless plasmas [5J. It 
has been shown that these particle-like excitations, which have topological character, 
interact not only with each other, but also with the curvature of the substrate [7]. 

In the case of ferromagnetic materials, cylindrical nanomagnets with a vortex as 
the magnetization groundstate have been considered as candidates to be used in logic 
memory, data storage, highly sensitive sensors [8J and cancer therapy [9]. Theoretical 
works have shown that the energy and stability of these topological excitations depend 
on the curvature of the ferromagnetic nanoparticle [10] and their dynamical properties 
are affected by both: the interaction with curve defects appearing during the fabrication 
of the nanomagnets and the presence of magnetic/non-magnetic impurities in these 
magnetic nanoparticles [T2] . 

In addition, vortices can appear like solutions to the continuous Heisenberg Model 
on two-dimensional systems. In Refs. [32] — [IS], the authors have used this model to 
analyse the dynamic and static properties of vortices and have shown that the energy of 
these excitations is closely linked to the characteristic length of the considered geometry. 
Besides, the vortex energy presents a divergence in simply-connected surfaces, which 
can be controlled by the development of an out-of-plane component in the vortex core, 
so called the vortex polarity. Soliton-like solutions have also been considered in the 
above cited works and it has been shown that their characteristic lengths depend on 
the length scale of the surface. For finite surfaces, fractional/half-soliton solutions have 
been obtained [151120]. 

Note that CMS in the presence of an external field exhibit changes in the behavior 
of collective modes and in their elastic properties. For example, an external magnetic 
field can be used to deform magnetoelastic metamaterials, which is achieved by pro- 
viding a mechanical degree of freedom so that the electromagnetic interaction in the 
metamaterial lattice is coupled to elastic interaction [2TJ. Furthermore, by combining 
the curvature effects with magnetic fields, the molecular alignment of flux lines of the 
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nematic director of nematic liquid crystals can be reoriented or switched between two 
stable configurations [22]. Lastly, it has been shown that the curvature of graphene 
bubbles can be controlled by applying an electric field [23] . 

In this context, Saxena et al. have considered the exchange and Zeeman terms in 
the magnetic energy calculations for cylindrical surfaces, and they have shown that the 
interaction of an external magnetic field with a cylindrical magnetoelastic membrane 
has a 2tt soliton-like solution, which induces a deformation (pinch) at the sector where 
the spins are pointing in the opposite direction to the magnetic field [23]. However, 
if a constant magnetic field is applied on an arbitrary curve surface, with rotational 
symmetry, the non-homogeneous double sine-Gordon equation (DSG) is obtained and 
numerical solutions are demanded [25J. However, from the coupling between the ex- 
ternal magnetic field and the curvature of a surface, analytical solutions are predicted 

In this paper, we show that the homogeneous DSG can be obtained for an ar- 
bitrary magnetically coated surface with rotational symmetry in the presence of an 
external magnetic field, provided the field is tuned with the curvature of the surface. 
The obtained solution consists in a 2ir soliton, which induces a deformation in the ro- 
tationally symmetric surfaces due the presence of the magnetic field. Furthermore, we 
apply the model to three specific surfaces: cylinder, catenoid and hyperboloid, in order 
to obtain the characteristic lengths of the 2ir soliton in these particular cases. The 
coupling between the magnetic field and the geometry of condensed matter systems is 
a interesting result and we believe that these results may guide future works in the 
control and manipulation of the shape and physical properties of magnetoelastic coated 
surfaces, superfluid helium, nematic liquid crystals, graphene bubbles, and topologi- 
cal insulators, once the results obtained here show that magnetoelastic surfaces may 
be deformed in a specific point by the application of a magnetic field varying from a 
maximum value, when p = 0, and tending to when p — > oo. 

This paper is outlined as follows: In section [21 we present the adopted model; the 
results and discussions are presented in Section El Section @] brings the application of 
the adopted model on three specific surfaces and in Section [5j we present the conclusions 
and prospects. 



3 



2 The model 



The total energy for a deformable, magnetoelastically coupled manifold is given by 
E = E mag + E e \ + E m _ e \, where E mSug , E e \ and E m . e \ are the magnetic, elastic and mag- 
netoelastic contributions to the total energy, respectively. The magnetic contribution 
is composed by the exchange, magnetostatic, anisotropy and Zeeman terms. In this 
work, we will focus our attention mainly on the exchange and Zeeman contributions 
to describe the magnetic properties of curved surfaces with rotational symmetry in the 
presence of an external magnetic field. In this case, the energy is well represented by 
the non-linear cr-model (NLcrM) on a surface in an external magnetic field: 



where m is the magnetization unit vector (m 2 = 1), dS is the surface element, B is the 
applied magnetic field, p is the magnetic moment, and g is the g factor of the electrons 
in the magnetic material. 

This model was previously considered for studying the properties of a circular 
cylinder surface in the presence of an constant axial magnetic field and the well known 
DSG [27] was derived from the Euler-Lagrange equations [21]. In that case, the authors 
have obtained a 2ir soliton like solution and have shown that a geometrical frustration 
appears, due to a second length scale introduced by the magnetic field. A surface 
deformation has been predicted at the sector where the spins point in the opposite 
direction to the magnetic field. On the other hand, unlike the cylinder case, whenever 
the above model has been considered on an arbitrary curved surface, it have yielded the 
non-homogeneous DSG, whose solution could be obtained only numerically [25]- Here, 
we show that the homogeneous DSG system appears for any rotationally symmetric 
surface, if the magnetic field is pointing in the z-axis direction and its strength is a 
function of p, that is tuned with the curvature of the substrate. That is, B = B(p)z, 
where p = p(z) is the radial distance from one point on the surface to the z axis. 

To proceed with our analysis, it will be useful to rewrite Eq.flTJ on a general ge- 
ometry with metric tensor gij, and rotational symmetry, described in cylindrical-like 
coordinate system. Thus, we have that g p ^ = g$ p = and the Hamiltonian (pQ) is 
rewritten as: 



where g pp and g^ are the contravariant metric elements. The magnetization unit vector 
is parametrized by m = (sin cos <£>, sin G sin $, cos G) and the rotationally symmetric 





H 




[gP" (d p Q) 2 + tf* sin 2 Q {d^f + gpB(p)(l - cos 9)] dpty, (2) 
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curved surfaces are parametrized in cylindrical coordinates system (p,<p,z). Besides, 
we have assumed cylindrical symmetry for the order parameter, that is, 0(p, 0) = 0(p) 
and $(p, <p) = $(0). It is important to note that, once the magnetic field is in the z- 
axis direction, the minimum energy for the Zeeman interaction is obtained when = 
and the maximum, when O = 7T. As it will be discussed after in this text, this fact 
forces the surface to deform (in order to minimize the total energy of the system) at the 
sector of the 2ir soliton where the spins point in the opposite direction to the magnetic 
field. 



3 Results and discussions 



We are interested in studying only static solutions of Heisenberg spins on curved 
surfaces in an external magnetic field. In this case, the Hamiltonian and Lagrangian 
of the system coincide, that is, H = L, where L = L[<d, $, G, $, t] is the Lagrangian of 
the CMS. Thus, in order to obtain the ground and excited states for this system, we 
must solve the Euler-Lagrange equations (ELE) for the Hamiltonian (j2J). Then, ELE 
yields: 



qpp 
—d 

g<M> p 



and 



nPP 

^ d P Q 

gH p 



sin 2 Qdd 



sin 26 + -^r-B'(p) sin 6 

2 y 1 



0. 




(3) 



(4) 



where B'(p) = gpB(p). 

Since we are considering surfaces with rotational symmetry, the parametric equa- 
tions associated with Eq.(3) and Eq.(4) can be written, in cylindrical-like coordinates, 
such as r = (p cos0, psin0, z(p)), where p is the radius of the surface at height z, and 
accounts for the azimuthal angle. In this case, the covariant metric elements are given 
by: 



and 



9 PP 



z' 2 + l, 



1 

gVV ' ' ' ■ gPP 

where z' = dz/dp. From the above equations, it is easy to see that 8$ C\ 
and Eq.fll]) is simplified to: 

sin 2 Qdl^ = 0. 



(5) 
= 0, 

(6) 



One can note that the above equation admits two kinds of solution: the first one 
is given by = nir, with (n = 0, 1, 2, ...). In this case, if n is even, the Zeeman term 



5 



in the Hamiltonian ( II ip is minimized, however, n odd yields a maximization of the 
magnetic field interaction energy, in such a way that these solutions must be unstable, 
since the spins would be pointing in the opposite direction to the magnetic field. Thus, 
the solution = (2n)-K is the ground state of Eq. (Tj[|). However, we are interested in 
the existence of excited states, so called 2ir solitons, which are a continuous transition 
between the two vacua, = and = 2tt, of this system. This class of solutions 
have been already obtained for the cylinder surface [21], however, analytical solutions 
have been also found for other geometries, provided the magnetic field is coupled with 
the curvature of the substrate [26J. Here, the simplest way to obtain topological 2ir 
soliton-like solutions is by considering the second simplest possible solution for the Eq. 
©, that is: 

$(0)=0 + o (7) 

where O is a constant of integration that does not influence the energy calculations. 
Then, we can rewrite the Eq.([3]) as: 

sin 20 

^6 = — + g M B'(p) sinQ, (8) 

where d£ = g^ / gppdp. 

From the analisys of Eq.®, there are two situations to be considered: firstly, when 
we take B(p) = 0, it yields the single sine-Gordon equation and a topological soliton, 
which satisfies the self-dual equations, is obtained [15J. This particular value for B(p) 
leads to the isotropic Heisenberg Hamiltonian, which has been previously studied on 
several curved geometries, e.g., torus [13], cylinder JT5J [16], sphere [17] . pseudosphere 
[18], cone [19], catenoid, hyperboloid [25]. It is important to note that, usually, the 
soliton presents a characteristic length scale (CLS) that depends on the geometrical 
properties of the underlying manifold. In general, the soliton's CLS appears in front 
of the sin(20) term, in Eq. flH]). However, this CLS is given by (y~g^~/ g pp ), which is 
embedded in the £ para meter. Thus, the soliton would have its CLS rescaled to one, 
e.g., if we take the cylinder case, the CLS of the soliton is equal to the constant radius 
p of this surface. However, when we are working on an infinite cylinder, the change 
of variable z —> z/p eliminates any dependence on p [T5] . 

The second case to be discussed is given when we consider an arbitrary magnetic 
field B(p) 7^ 0. In this case, if the magnetic field is constant, and an arbitrary curved 
surface with rotational symmetry is considered, Eq.flS} leads to the non- homogeneous 
DSG and numerical solutions are demanded [25]. However, it can be noted that if 
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the external field is coupled with the surface curvature, in the form B'(p) = g^B' , 
where B' = gpB(p ) and p is the surface radius at z = plane, Eq.© yields the 
homogeneous DSG: 



^ = ^ + 5 ; S i„e. (9) 



sin 26 

y 

Assuming that 0(— oo) = and 6(+oo) = 2n, the solution for the above eqquation 
can be given by: 

e( « = 2tm " 1 (^)' (10) 

where p 2 = 1/B' Q and C = P B /{1 + P 2 ) 1/2 - Then, one can conclude that the homoge- 
neous DSG must be obtained for the Heisenberg spins on curved surfaces in an external 
magnetic field only if the magnetic field is tuned to the surface curvature, varying in 
function of l/p(z). 

Eq. (|T0"|) represents a 2-7T soliton, which is a topological excitation belonging to the 
second class of the second homotopy group, whose CLS is given by (. Note that the 
increasing of B' leads to the decreasing of the CLS of the soliton causing it to remain 
confined in smaller regions of the surface. It is also easy to see that ( — > when 
B' Q — > oo. Note also that, if we put B' = in Eq. flTDl) . a it soliton is not obtained, 
as it happens when we consider Eq.Q with B' = 0. This is explained because these 
two solutions belong to two different homotopy classes and cannot be transformed one 
to another by continuous transformation or by limiting process B' — > 0. It can also 
be noted that, by introducing a new CLS in the system, the magnetic field induces 
a geometrical frustration, and the soliton, which had its CLS rescaled to one, must 
choose a new CLS, given by (, that is smaller than the length p B , introduced by the 
magnetic field. 

In order to calculate the energy of the 2ir soliton, we can rewrite the Hamiltonian 
(ED as: 



H' = H 1 + H 2 = 2vr 



/OO POO 1 

[(<9 ? 6) 2 + sin 2 6] d£ + / — (1-cosG) 
-OO J —OO Pa 



-co r B 

where Hi is the part of the Hamiltonian that does not have dependence on the magnetic 
field and H 2 is the Zeeman term. In this case, the energy of the soliton is evaluated to 
give: 

7 i\ 1/2 i 

E a = 8vr 1 + — + — sinh" 1 p B , (12) 

V p\) pI 



that is larger than the minimum energy for the homotopy class with winding number 
Q = 2, that is to say, E 2iv s = 8n. Only in the limit p B —> oo, we get E s — > 811 and 
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£ — > 1. Once p B = y/l/B'(p ), the energy of the soliton is associated with the strength 
of the magnetic field at the z = plane. In order to release the geometrical frustration 
introduced by the magnetic field, an elastic surface will deform to decrease the radius 
in the region where the soliton is centered. This will minimize the energy of the two 
terms that do not depend on the magnetic field in the Hamiltonian (|lip . Furthermore, 
by considering only the first term the development of ELE leads to: 

fl?e=*£° (is) 

as expected. This is the single sine-Gordon equation, whose solution is: 

6 = 2 arctan (e c ) , (14) 

which represents a 7r soliton, interpolating the two minima, = and = tt, and 
which belongs to the first class of the second homotopy group. As it was mentioned 
before, the 7r soliton solution can not be obtained if we take the limit B — > in 
the solution (fTUl) . as well as we do not obtain the tc soliton energy, predicted by the 
Bogomolnyi inequality [28] by taking B = in the Eq. (I12j) . 

Now, we will focus on the last term in the Hamiltonian H', which involve the 
external magnetic field: 

,oo 1 

H 2 = 2n — (l-cose)df. (15) 
J-oo P 2 B 

For the solution (flQl) . H 2 will take the form: 

/OO -I 1 /»oo 

-ooP B [1 + (l-C 2 )smh 2 |J J-oo 

This part of the Hamiltonian will decrease in energy if we locally decrease £, but keep 
it constant in Eq.f llOp . Therefore, for a magnetic field tuned with the surface, if we 
keep the surface cross section at £ ± oo, the surface will deform in the region of the 
soliton. In the case of small magnetic field, we get: 



E s = 8tt 



l + -^(l + 21n2p £ 



D 



(17) 



and for large magnetic field, we get E s = 8n/p B . 

As expected, the highest contribution for the magnetic energy density due to the 
interaction with the external magnetic field comes from the sector where the spins 
point in the opposite direction to the field. In this way, the 2tt soliton would like to 
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collapse and eliminate the region where the spins are opposite to the magnetic field. 
However, to align the spins, the system would fall in a tt soliton sector of the second 
homotopy group, which has infinite energy due to the interaction of the spins at +00 
with the magnetic field. Thus, there is a hard-core repulsion between the two tt solitons 
in the system by virtue of the curvature. Due to the introduction, by the magnetic 
field, of a new CLS into the system, a geometric frustration appears in the problem 
and the energy in Eq. (fl2l) is grater than 871". Then, the surface will try to deform 
in order to decrease its magnetic energy until it reaches the second homotopy class 
[21]. Obviously this deformation will cost elastic energy to the system, so the real 
deformation will balance between the gain in magnetic energy and the loss in elastic 
energy.This result suggests that we could, at first sight, deform an elastic membrane, 
coated with magnetic material, by using a variable magnetic field, whose strength 
would decrease with the increasing of the surface's rotational radius, p. 

In the absence of an external field (B = 0), the self dual equation for the Eq.Q 
is d^Q = ±sinB. For the DSG, we can write: 



The function in square brackets is greater than one for any value of 0(£) and the self 
dual equation is not satisfied. Therefore, in order to satisfy the self dual equation and 
to lower the energy, the surface would deform. Note that all solutions of the Eq. (1181) 
satisfy the Euler-Lagrange equation, but not vice versa. The 2tt soliton lattice of the 
double sine-Gordon equation in two regimes, p B < 1 and p B > 1, are similar to that 
found on the cylinder [21]. 




4 Some specific cases 

27r-soliton solutions should be found on infinite surfaces with rotational symmetry 
and the CLS is associated with the geometry of the substrate. In this case, we will 
apply our results on some infinite surfaces with this kind of symmetry in order to find 
the soliton CLS associated to each one of them. Here, we will consider three particular 
geometries: the cylinder, the catenoid and the hyperboloid, whose shapes are shown in 
Fig. HI It can be noted that in spite of the shape similarity between the catenoid and 
hyperboloid surfaces, their geometric properties are different, as it will be discussed 
soon after. It is also important to note that the 27r-soliton, given by the solution 



<%0 = ±sin© 1 + 



4 sin 2 (9/2) 
p 2 sin 2 O 

' B 



(18) 
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Figure 1: [Color online] From left to right, one can see the surfaces of cylinder, catenoid 
and hyperboloid. These are curved surfaces with rotational symmetry. In order to 
obtain a 27r-soliton on these surfaces, the magnetic field must be coupled with the 
curvature. In the case of the cylinder, the magnetic field is constant, while that one 
for the catenoid and hyperboloid varies with p, having maximum value in the z = 
plane, and vanishing its strength in p — > ±00. 

ffTDl does not appear on finite surfaces, e.g., sphere, torus and/or barrels. In these 
cases, fractional solitons must be found, because the spin sphere will not be completely 
covered two times, as it should be in order to get an excitation belonging to the second 
class of the second homotopy group. In this way, topological arguments can not be 
used to ensure the stability of this excitation. 

Thus, the first case to be analysed is the surface of the cylinder, which was pre- 
viously studied by Saxena et al [21] • In that work, the authors have predicted the 
appearing of a 27r-soliton on the cylinder in the presence of a constant magnetic field 
and have found its CLS as a function of z. The cylinder presents null Gaussian curva- 
ture and its mean curvature is given by M = l/2r. In this case, we have that z(p) = z 
is a constant, thus, the cylindrical coordinates parametrization leads to the metric 
elements: 

9t = 9 2 = r 2 and g% = 1, (19) 

where r is the radius of the cylinder. It is easy to see that the magnetic field to be 
applied on the surface is B'(p) ci = (l/r 2 )B' , which is constant. The CLS is easily 
calculated to be: 

£ cyl = lnr. (20) 

As expected, the result here obtained would be the same as that given in Eq.(4) of Ref. 
[21] if we had considered z instead p as the variable of the problem. 
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From the analysis of Eq. OH]), we could conclude that the homogeneous DSG can 
be obtained from the Heisenberg model in an external field for any other rotationally 
symmetric surface, besides the cylinder, if the magnetic field has its strength varying 
with the metric element g^. On the other hand, if we apply a constant magnetic field 
on an arbitrary curved magnetoelastic surface, the solutions to the model described by 
the equation ([1]) lead to the inhomogeneous DSG, which must be solved numerically 
[25] . We continue our analysis, considering two surfaces with negative varying Gaussian 
curvature, namely the catenoid and the hyperboloid. 

The catenoid is a minimal surface, with its mean curvature equals zero everywhere. 
However, unlike the cylinder surface, its Gaussian curvature does not vanish. Indeed, 
the the catenoid presents a negative Gaussian curvature, which is given by: 

* = (21) 

Another difference of the catenoid in relation to the cylinder is the act that the radius 
of the catenoid is not constant, but it varies as a function of z, as follows: 

p(z) = rcosh(z/r). (22) 

The catenoid can be parametrized by (p cos (f), p sin <fi, r cosh -1 (p/r)), where r is the 
radius of the surface at the z = plane. This parametrization leads to: 

■?_ 
— r 

Then, in order to obtain the homogeneous DSG, the magnetic field to be applied on 
this surface, will be -B'(p) cat = (l/p 2 )B' . Once p is given by Eq. ff22l . the magnetic 
field is not constant, but it varies in function of z. Indeed, for p — > oo, the field tends 
to zero and its maximum value (B'(p = r) cat = (l/r 2 )B' ) must occur at the point 
p(0) = r. Finally, it is easy to calculate the CLS on the catenoid in an external varying 
magnetic field, which admits a 27r-soliton solution with: 

£ cat =ln[2(p +v /7^)l. (24) 



'^P 2 and C = ~T-rr ( 2 3) 



The last surface on which our model will be applied is the hyperboloid, which has 
a shape similar to the catenoid, however, while the catenoid has mean curvature zero 
everywhere, the hyperboloid has variable mean and Gaussian curvatures [21] • A one- 
sheeted hyperboloid can be parametrized by (p cos <p, p sin 0, (b/r) \J p 2 — r 2 ) , where, 
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again, r is the radius of the surface at z = and b is a multiplicative parameter that 
accounts for the height of the surface. This parametrization leads to: 



44 -P and # pp - . [Zb) 



p 2 {b 2 + r 2 ) - r 4 

From now on, in order to obtain analytical solutions, we will use the biharmonic coor- 
dinate system (BC), in which b = r, to describe a particular kind of hyperboloid [30] . 
here called polar hyperboloid. In this case, the metric element simplifies to: 

?n 2 - r 2 

€ yP = jr—^ (26) 

and, to obtain the homogeneous DSG, the magnetic field must be given by B'(p) h = 
(l/p 2 )B' . Thus, as well as in the catenoid case, in order to obtain a 27r-soliton, the 
magnetic field must be p dependent, having its maximum value, given by B\p = 
r )h yp = (l/ r2 )-^cb & t the z = plane. It is easy to show that the magnetic field applied 
on the hyperboloid described by BC varies with p = \/r 2 + z 2 , then, B'(p) hyp — > 
when z(p) — > ±oo. Finally, the polar hyperboloid also admits the solution given by the 
equation ( TlOl) . however, the £ parameter for this surface is given by a large and tedious 
expression which we will not consider here. Despite the fact that lim p _ >00 B(p) = 0, 
it is important to note that the total magnetic flux through these infinite surfaces di- 
verges as In p at infinity for the catenoid and hyperboloid, but remains constant and 
depending on r for the cylinder. 



5 Conclusions and prospects 

In conclusion, we have shown that the Euler-Lagrange equations derived from the 
continuum approach for classical Heisenberg spins on a surface with rotational symme- 
try in an external magnetic field, pointing in the z axis direction, is the homogeneous 
DSG provided the field is coupled with the curvature of the surface. We have found a 
single 27r-soliton like solution for this model for an arbitrary surface. Surface deforma- 
tions at the sector where the spins point in the opposite direction to the magnetic field 
were predicted, in order to diminish the energy, which is greater than that predicted 
by the Bogomolnyi's inequality. 

We have applied this model on three specific surfaces: cylinder, catenoid and hy- 
perboloid, which have different geometrical properties, characteristic length scales and 
Gaussian curvatures. As expected, each considered surface admits a 27r-soliton, but the 
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magnetic field which gives the homogeneous DSG presents different behavior on each 
one (in each case B is a function of p{z)). While the field to be applied on the cylinder 
must be constant, the field which leads to the homogeneous DSG on the catenoid and 
the hyperboloid is inversely proportional to p(z). Since the characteristic length of 
the soliton depends on the strength of the field, the magnitude of the deformation can 
be controlled by the magnitude of the applied magnetic field. Finite surfaces admit 
fractional 27r-solitons, which are not topologically stable. 

Our results can play an important role for experimental studies where magnetoe- 
lastic effects here predicted could be observed. For example, a membrane coated by 
a magnetic material could be deformed by using this technique. Additional theoret- 
ical work which also includes surface tension could clarify the deformation dynamics 
of magnetic coated membranes and/or magnetoelastic metamaterials pljj. Finally, we 
believe that this issue may be relevant for future studies on manipulation and control 
of the shape and physical properties of nematic liquid crystals, curved graphene sheets 
and topological insulators. 
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